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I. INTRODUCTION 



The theoretical tools for the description of exclusive processes are the hard-scattering am- 
plitude which describes the pro cess- dependent quark-gluon interaction within perturbative 
QCD and the probability amplitude for finding the lowest twist quark state in each hadron. 
The total amplitude is represented by the convolution of these two parts , assuming factor- 
ization of highly off-shell or large transverse momentum regions of phase space from regions 
of low momenta necessary to form bound states. Higher-twist components, corresponding to 
a higher number of partons (quark-pairs and gluons) are supposed to be suppressed by pow- 
O . ers of the momentum transfer Q^. Recent progress 0,^] in Sudakov-suppression techniques 
^ ! provides support for the conjectured infrared protection of the perturbative picture. The 
^ I focus in this talk will be on recent theoretical developments on exclusive reactions involving 
P^l nucleons. 

D , 

^ . II. PERTURBATIVE ASPECTS 

■ "^^^ Q^-dependence of hadronic wave functions is determined by the renormalization 

^ ' group equation. To leading order, the dynamical evolution of the lowest twist hadronic 
distribution amplitude, which is the hadronic wave function integrated over transverse mo- 
menta up to a resolution scale Q^, is described by the one-gluon exchange kernel subsumming 
leading logarithms of ladder graphs. Specifically, the nucleon distribution amplitude $n is 
obtained as the solution to the integrodifferential equation 

{q'^, + ^]hxuQ') = ^ l\dy]V{xum)Hm,Q') (i) 

[Jg [fix] = dxi dx2 dxs6{l — Xi — X2 — X3)], where Cp and Cb are the Casimir 
operators of the fundamental and adjoint representations of SU{3), respectively, and jS is 
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the Gell-Mann and Low function. The leading-order expression for the integral kernel V has 
been calculated in 

To solve the evolution equation (1), ^nI^^i, Q^) has to be expanded in terms of the eigen- 
functions of the integral kernel V. These eigenf unctions correspond to three-quark operators 
which are multiplicatively renormalizable, i.e., to operators with definite anomalous dimen- 
sions [|]. The latter can be determined by diagonalizing the evolution equation within some 
appropriate basis. To this end, one expresses the solution of (1) in the form 

<i>N(xi,g2) = <^^^:,^)f2BMxi){^f^y\ (2) 



n=0 



where {$n}=o are orthonormalized eigenfunctions of the interaction kernel within a truncated 
basis of Appell polynomials of maximum degree M. ^a.s{xi) = 120x1X2X3 is the asymptotic 
form of the nucleon distribution amplitude Because the 7n are positive fractional numbers 
increasing with n, higher terms in (2) are gradually suppressed. The Appell polynomials 
are polynomials of two independent variables, say xi and X3. Thus one can expand {$n} in 
terms of the polynomial basis {x\x3}?j'^Q: 



l>n = 5^a:;xixi. (3) 

ij=0 



Then defining moments 



^N^f^') = f[dx] x[ x° xi $n(x,,/x^), (4) 
Jo 

the expansion coefficients -Bn can be formally determined by inverting (2) to obtain 

An explicit calculation of the expansion coefficients -B„(/x^) involves the orthonormalization 
of polynomials with two variables — a problem with no unique solution because the procedure 
depends on the order in which it is performed. To determine the anomalous dimensions 

7n = (^1^ + 2^711^ j 5 one has to compute first the zeros t]^ in the characteristic polynomial 

that diagonalizes the evolution kernel. Such a program has been carried out in and a 
complete eigenfunction basis has been constructed up to polynomial order 9. All previous 
calculations had been restricted to maximum order 3 . It is noteworthy that up to order 
7, we have performed the diagonalization of the evolution kernel analytically. Below order 
3, our anomalous dimensions coincide with those computed by Peskin 0; those of order 3 
confirm the recently published (numerical) estimates of [^. Our results are listed in Table I; 
for more details we refer to 0. The large number of computed eigenvalues — a total of 54 
corresponding to 29 symmetric and 25 antisymmetric eigenfunctions — enables the evaluation 
of a well-defined pattern (Fig. 1). The trend line of this pattern follows the empirical power 
law (solid line) 7n = 0.37 0(n)°'^^^. There is, certainly, no evidence that the ensuing global 
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behavior of baryon anomalous dimensions coincides with that of mesons (dotted hnes), as 
prematurely suggested in [^. The proposed method can be used to consistently generalize 
the ansatz for $n to higher orders by retaining the correct evolution behavior 

The calculation of moments of the nucleon distribution amplitude resides on nonpertur- 
bative techniques and will be considered in the next chapter. 
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TABLE I. Orthogonal eigenfunctioiis $„(.ti, 2:2, .T3) = u,^;x;^j.3 
equation (represented by the coefficient matrix d'^f with aj!^ = S'n a^.; 
normahzation is given by Jq [dx] xia;2X3 = (A^n)~^ (^fcn- 



of the nucleon evolution 
a^2 = for ah n). The 
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FIG. 1. Eigenvalues of the evolution equation for symmetric (left-hand side) and antisymmetric 
eigenfunctions (right-hand side). 



III. NONPERTURBATIVE ASPECTS 



The derivation of the nucleon distribution amplitude from QCD is intimately connected 
with confinement and employs nonperturbative methods, e.g., QCD sum rules |Q, lattice 
gauge theory ^ or the direct diagonalization of the light-cone Hamiltonian within a dis- 
cretized light-cone setup [|10[- To determine the moments ^(^^"'^"-a)^ ^ short-distance operator 



product expansion is performed at some spacelike momentum /i^ where quark-hadron duality 
is valid M. One considers {z is a lightlike auxiliary vector with = 0) 



^z.-^) Mz,.p) 



Zi=0 i=l 







(z-p) 



(6) 



and evaluates correlators of the form [pj] , pr^ 



/("i"2"3''")(g,2) =z / d^xe^^-^ < f]|T(O('^i"2"^)(0)O^r^(a;))|f^ > (2-7) 



/77' 



(7) 



where the factor 
and 



7)^^, serves to project out the leading twist structure in the correlator. 



(8) 



are appropriate three-quark operators containing derivatives. Their matrix elements 

< fi|o;^"^''2"«)(o)|P(p) >= /jv(2-p)"^+"'+"^+^A/'^o("^"2"«) (9) 



are related to moments of the covariant distribution amplitudes ||13| V , A, and T: $N(a^i) = 
V{x,) - v4(xi), $n(1,3,2) + <I>n(2,3,1) = 2T(1,2,3) with V(1,2,3)=V(2,1,3), A(l,2,3)=- 
A(2,l,3), and T(1,2,3)=T(2,1,3). Here /n denotes the "proton decay constant". 

On the basis of such QCD sum-rule calculations, useful theoretical constraints on 
the moments of nucleon distribution amplitudes have been obtained and various mod- 
els |18| JTT|JI^JI^JI5[| have been proposed. Examples of physical observables calculated from 



these models are given in [^,|rT|,|I^] and more recently in |T3| , P!7| -[TD[] . Our project differs from 
previous ones in that we use a "hierarchical" x^-criterion to parametrize the deviations from 
the sum-rule constraints ||2^ . This affords to determine optimized versions of previous model 
distribution amplitudes as well as to find a new solution (we labeled heterotic) which 



hybridizes morphological and dynamical features of COZ-type |Tl[] and GS-type ampli- 
tudes providing results corroborated by the available data. The "hierarchical" treatment of 
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the sum rules takes into account the higher stabihty of the lower-level moments [|T6[ and does 
not overestimate the significance of the still unverified constraints |11| for the third-order 



moments. The simple but important assumption is that the model space can be safely trun- 
cated at states with bilinear correlations of fractional momenta because adding higher-order 
contributions should only refine the initial approximation. The advantage of our method 
becomes apparent by taking a more global approach, i.e., looking for solutions on the scale 
of the whole validity range of the sum rules. This treatmant leads to a characteristic series 
of local minima shown in Fig. 2. 



FIG. 2. Large scale pattern of nucleon distribution amplitudes complying with existing QCD 
sum rules. 

They constitute a pattern characterized by a smooth and finite orbit which is com- 
pletely specified by joint values of -B4 and R = IG^I/G^. This striking scaling behavior 
seems to pertain to solutions p2|,|7[| which incorporate higher-order eigenfunctions (Appell 
polynomials) — not used in our fit procedure (see inset in Fig. 2). The robustness of the 
fiducial orbit suggests that many of the features of nucleon distribution amplitudes that 
seemed unrelated actually fit together into a coherent overall structure, hence leveraging 
our knowledge of specific characteristics into a more general context. Isolated samples in 
the (i?4, R) plane are relegated to spurious solutions, either because they exhibit unrealistic 
large oscillations in the longitudinal momentum fractions [^] (stars) or because they yield a 
wrong evolution behavior for the nucleon form factors (light upside-down triangles). The 
presented curves are fits to the local minima of the COZ sum rules (+ labels) and to a com- 
bined set of KS and COZ sum rules (o labels) restricted within the intervals 0.104 0.4881 
and 0.0675 ^ 0.482, respectively: R = 0.437338 - 0.00601654 - 0.000176S| (solid line) and 
R = 0.431303 -0.00752^4 -0.00024151 + 3.851221 x lO-*^^! (dotted line). Cross-type solu- 
tions show up across the fiducial orbit and they seem to undergo a complete metamorphosis 
as they pass from COZ-type solutions {R =< 0.5) to the heterotic one (smallest possible 
value of R still compatible with the sum-rule constraints). The passage between these two 
types of solutions is smooth and the changed shapes follow an orderly sequence of grada- 
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tions characterized by comparable values of x^- AH solutions are organized by the same kind 
of dependence between R and -B4. The solution denoted Het' , past the COZ-cluster, has 
"mirror" -image characteristics relative to the heterotic solution pll|18| but is unstable. 



IV. SUMMARY AND CONCLUSIONS 

We have attempted to provide a unified view of nucleon distribution amplitudes derived 
on the basis of QCD sum rules. The pattern that emerges is not one of widely dispersed 
solutions across the parameter space, but one of a smooth and finite curve (an orbit) in the 
(i?4, R) plane, as underlined in Figure 2. On the perturbative side, a systematic theoretical 
procedure to calculate orthonormalized eigenfunctions of the leading-order nucleon evolution 
kernel has been developed and a complete set of the first 54 terms has been explicitly 
evaluated. 

This work was supported in part by the Deutsche Forschungsgemeinschaft. 
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